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Spontaneous stratification of granular mixtures has been reported by Makse et al. [Nature 386, 
379 (1997)] wlien a mixture of grains differing in size and sliape is poured in a quasi- two-dimensional 
heap. We study this phenomenon using two different approaches. First, we introduce a cellular 
automaton model that illustrates clearly the physical mechanism; the model displays stratification 
whenever the large grains are rougher than the small grains, in agreement with the experiments. 
Moreover, the dynamics are close to those of the experiments, where the layers are built through a 
"kink" at which the rolling grains are stopped. Second, we develop a continuum approach, based on 
a recently introduced set of coupled equations for surface fiows of granular mixtures that allows us 
to make quantitative predictions for relevant quantities. This approach includes amplification (i.e., 
static grains entering the flow of rolling grains) , a phenomenon neglected in the cellular automaton 
model. We study the continuum model in two limit regimes: the large flux or thick flow regime, 
where the percolation effect (i.e., segregation of the rolling grains in the flow) is important, and the 
small flux or thin flow regime, where all the rolling grains are in contact with the surface of the 
sandpile. (1) In the thick flow regime, where most experiments are carried out, the flowing grains 
are segregated in the rolling phase; as they are flowing down, the large rolling grains are convected 
to the top of the layer, and only the small rolling grains interact with the sandpile. We include 
this effect in the continuum model and flnd results very close to the experiments. (2) In the thin 
flow regime, we find interesting results that are close to the thick flow limit. However, due to the 
presence of cross-amplification processes, we flnd a small regime that shows stratiflcation when the 
small grains are slightly rougher than the large grains, but stratification is much more pronounced 
if the large grains are rougher. We study in detail the dynamical process for stratification, where 
the layers are built through a "kink" mechanism, and find the dependence of the size of the layers 
on the parameters of the system. We flnd that the wavelength of the layers behaves linearly with 
the flux of grains. We also flnd a crossover behavior of the wavelength at the transition from the 
thin flow to the thick flow regime. We obtain analytical predictions for the shape of the kink giving 
rise to stratiflcation as well as the proflle of the rolling and static species when segregation of the 
species is observed. Finally, we compare our approach to previous ones. 



X 



I. INTRODUCTION 

Segregation is a prominent example of the peculiar 
properties of granular matter ||l|-^. For example, shak- 
ing a container filled with two types of grains of different 
sizes leads not to mixing — as in liquids — but to segrega- 
tion, with the large grains on the top of the container 
and the small grains on the bottom. This striking be- 
havior and the importance of mixing problems from a 
technological point of view have led to a broad interest 
in granular materials in the physics community. 

Several types of segregation have been investigated, 
namely, segregation by vibrating a container (the "Brazil 
nut effect") [p^-p4[, and segregation in rotating cylin- 
ders, where the segregation occurs through surface flow 

^5HT1. 



Segregation can also been obtained in the absence of 
any periodic oscillation by simply pouring a mixture of 
grains of different sizes onto a pile. One experimental set- 
up, which has attracted much recent attention, consists 
of a quasi-two dimensional cell or vertical Hele-Shaw cell 
where a mixture of grains is constantly poured next to 
one end of the cell [Fig. |^(a)] . When a mixture of small 
and large grains is poured into the cell, a pile builds up 
and the small grains are observed to segregate near the 
top of the pile and the large ones near the bottom of the 
pile |p0|-p6t . This segregation is due to the different grain 
sizes, because large grains roll down more easily on top 
of small grains than small grains on top of large grains. 



It has been recently observed ]27|-|30[| that if the shapes 
of the species are sufficiently different, a novel type of 
segregation can take place: when such a mixture of small 
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rounded grains and large rough grains is poured between 
two vertical slabs, a spontaneous stratification of the mix- 
ture in alternating layers of small and large grains par- 
allel to the pile surface is observed. Additionally, there 
is an overall tendency for the large and small grains to 
segregate into different regions of the cell. Spontaneous 
stratification is also found for more than two species . 
The layers are then ordered forming a sequence parallel 
to the surface pile (from bottom to top, small-medium- 
large, small-medium-large for 3 species). 

This phenomenon could possibly be relevant in differ- 
ent fields, such as geomorphology |^,^. For example, 
stones coming from sand dune solidification (sandstone) 
present successive alternation of layers of different types 
of grains. This regularity cannot be explained by peri- 
odic sedimentation. But one could imagine that the sand 
dune was built with sand brought by the wind and flow- 
ing down the slip- face of the dune Q| , where stratification 
could appear as in the experiments of [p7| . This phe- 
nomenon may also have important consequences in in- 
dustry where processing and transport of initially mixed 
grains could lead to stratification. 

Spontaneous stratification has been observed in p7| af- 
ter a transient regime and only when the large grains are 
rougher than the small grains. When the large grains are 
the smoothest stratification does not occur but only the 
segregation of the mixture in different regions of the cell 
occurs, with the large and rounded being at the bottom of 
the cell, and the small and rough at the top. These phe- 
nomena can be qualitatively understood as the growth 
of an instability due to the competition of two opposite 
effects, (i) On the one hand, if the grains have the same 
shape (spherical for example) or, more precisely, the same 
repose angle, then the large grains will roll more easily 
on the top of the small ones than the reverse. This will 
lead to segregation, with the small grains at the top of 
the pile and the large at the bottom, (ii) On the other 
hand, if the grains have the same size but different shape, 
the rougher grains will have a larger friction coefficient 
than the other grains, and we expect segregation with the 
rougher grains at the top of the sandpile. In the case of 
large grains rougher than small grains, those two effects 
will compete, giving rise to an instability, and then to 
spontaneous stratification. However, in the case where 
the smaller grains are rougher than the large grains, the 
two effects will contribute in the same way, and the seg- 
regation pattern will be stable, without the occurrence 
of stratification |Q . 

In this paper, we study granular flow in the spirit of 
P^-^, where one assumes a sharp distinction between 
the rolling or fluid phase composed of the layer of grains 
flowing down and the static phase or bulk composed of 
grains forming the sandpile. The essential feature to un- 
derstand sandpile formation is then to describe how the 
rolling grains are flowing down and how they interact 
with the static grains. In this context, it is important 
to distinguish two different regimes, the thin flow regime 
and the thick flow regime, where the interaction between 



rolling grains and static grains are very different. More- 
over, experiments have been done in both regimes, so 
that both regimes are important to study. 

Here we propose a mechanism to understand sponta- 
neous stratification and segregation in two-dimensional 
silos. We adopt two different approaches. First we intro- 
duce a cellular automaton model to give insight into the 
proposed mechanism. This model displays stratification 
as soon as the large grains are rougher than the small 
ones, in agreement with experiments pTj . The dynam- 
ics are also very close to the experiments , where the 
layers are built through a "kink" going uphill where the 
rolling grains are stopped. 

Second, we study a continuum model which is built 
on a phenomenological formalism describing granular 
flow on sandpiles. This formalism was introduced by 
Bouchaud et al. jsj] for the case of a single species, and 
recently generalized by Boutreux and de Gennes (BdG) 
p6| to the case of a mixture of two species. The essen- 
tial feature to understand sandpile formation is to de- 
scribe the interaction between the rolling or fluid phase 
and the bulk, as well as any segregation which may oc- 
cur in the rolling phase itself. In contrast to the cellular 
automaton model, this formalism includes amplification 
processes, i.e., the conversion of static grains at the sur- 
face of the sandpile into rolling grains by the flow (this 
feature is essential to understand avalanching when one 
tilts a sandpile above the maximum angle of repose, or 
to understand avalanche dynamics in rotating drums). 

In the large flux regime, the rolling grains form a thick 
phase and the grains are kinematically segregated in the 
rolling phase, an effect called kinematic sieving, free sur- 
face segregation, or percolation [Fig. |^(b)] ]2^ , ^ , p8[ . 
Due to this phenomenon the large grains in the rolling 
phase are found to rise to the top of the rolling phase 
while the small grains sink downward through the gaps 
left by the motion of larger grains in the rolling phase. 
Thus, small rolling grains form a sub-layer underneath 
the sub-layer of large rolling grains. Then only the lower 
grains of this layer interact with the sandpile. Thus, a 
segregation effect that we refer to as percolation takes 
place inside the rolling phase where the large grains are 
convected along the flow to the top of the layer. Thus, 
if small rolling grains are present, they will preferentially 
be on the bottom of the layer, preventing the large grains 
from interacting with the sandpile. 

We will develop a continuum formalism that includes 
the percolation effect, and show that stratification arises 
in a similar way as in the experiments. When the large 
grains are rougher, stratification is made of layers spread- 
ing all along the sandpile; the layers of the two types of 
grains being of the same size for an equal volume mixture 
of grains. When the small grains are rougher than the 
large grains, we observe the complete segregation of the 
mixture but not stratification, with the small and rough- 
est grains being found at the top of the pile in agreement 
with the experiments. A thick rolling phase is a condition 
for the percolation effect to take place. However a large 
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difference in size is also necessary to observe the percola- 
tion effect: we expect d2/d,i ^ 1.5, where di and d2 are 
the typical size of the small and large grains, respectively. 

In the small flow regime each rolling grain is always 
in contact with the sandpile and so interacts with it. In 
this case a thin rolling phase is expected — of the order 
« d2 — and percolation effects are expected to be irrel- 
evant. However strong segregation effects are expected, 
as in the case of percolation, since we will consider the 
case of large size difference between the grains. The 
large grains are not being captured on top of static small 
grains (although they interact with them) because of the 
large difference in size, while the small grains are eas- 
ily captured on top of the large static grains. Using the 
theoretical formalism adapted to this case, we will show 
that stratification is observed when the large grains are 
rougher than the small grains, in agreement with the ex- 
periments. 

Cross-amplification processes — whose influence was 
greatly reduced in the case of thick flows due to percola- 
tion effects — now plays an important role since all grains 
are interacting with the pile and cross-interactions of the 
type small/large grains are expected to occur at the fluid- 
bulk interface. Due to cross-amplification, the transition 
from stratification to segregation does not occur sharply 
when the small grains become the roughest. We find a 
regime where the small grains are slightly rougher than 
the large grains, which also gives rise to stratification. 
When the small grains are more rounded than the large 
grains, then we find the segregation of the grains in agree- 
ment with experiments. 

In general, we find that the thick flow regime and the 
thin flow regime show similar results regarding the strat- 
ification and segregation of the species. This is valid pro- 
vided that the grains differ appreciably in size (in term 

of the size ratio we expect ^2/^1 ~ 1.5 |^,^). Then 
we study analytically a mechanism for stratification and 
segregation valid for both regimes. We study the "kink" 
mechanism, and we obtain predictions for the wavelength 
of the layers as a function of the different parameters of 
the system. In particular we find that the wavelength 
of the layers is proportional to the flux of adding grains 
(i.e., proportional to the thickness of the rolling phase), 
and we also find a crossover behavior at the thin-thick 
flow transition. Finally we discuss our results in regard 
to previous approaches to segregatio n [p6| and stratifica- 
tion and to the experiments p7|~[30|| . A short report 
of some of the results presented here has been published 
in 13^1. 



models have been used before for modeling the complex 
behavior of granular materials. The role played by the 
angle of repose has been used to understand the complex 
dynamics of granular flows of mixtures [p|, ^4|j3^j39[| . The 
concept of angle of repose normally involves a macro- 
scopic definition — i.e., the angle of repose is the angle of 
the pile surface after the grains are poured onto a heap. 
Here, we will define a microscopic angle of repose as the 
maximum angle at which a rolling grain is trapped at the 
surface of the pile. This definition must be understood 
in terms of a coarse grained distance along the surface of 
the pile — of the order of a few grain sizes — over which all 
the quantities involved in the model are averaged. 
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FIG. 1. (a) Typical experimental set-up to study sponta- 
neous stratification and segregation in quasi-two-dimensional 
cells, along with the different quantities defined in the text, 
(b) Schematic representation of the percolation effect in the 
rolling phase. We pour continuously a mixture of large and 
small grain in a two-dimensional cell and observe a steady 
thick flow of grains. The percolation effect consists of the size 
segregation of grains in the rolling phase: large grains are ob- 
served to rise to the top of the rolling phase, and small grains 
drift to the bottom of the rolling phase. Due to the percola- 
tion effect, the small grains are the flrst to be captured at the 
pile surface, resulting in the segregation of the mixture in the 
bulk. 



II. CELLULAR AUTOMATON MODEL 

A. Motivation and Definitions 

We first introduce a simplified discrete model to under- 
stand the mechanism leading to stratification. Discrete 



In the case of a single-species sandpile, as sand is added 
to the sandpile we consider a critical angle of repose Or- 
When the local angle of the pile is smaller than 9r the 
rolling grain is stable and there is no flow. When the 
local angle of the pile is larger than Or the rolling grain 
is unstable and flow occurs. 
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We develop a simple discrete model, based on the idea 
that grains flow when the local angle of the sandpile is 
larger than the angle of repose. We consider that the 
angle of repose depends on the type of rolling grains and 
also on the composition of the surface. Thus, in the case 
of granular mixtures of two different species, we consider 
the existence of four different generalized angles of repose. 

The sandpile is built on a lattice, where the grains have 
the same horizontal and vertical size as the lattice spac- 
ing. The two species will be distinguished according to 
different generalized angles of repose. In general, we will 
call species type 1 to the small grains, and species type 
2 to the large grains. The species can also have rough 
or smooth surface properties, which, together with the 
different size, will define the different angles of repose of 
the species. Following ||3^-|36t, we regard each grain as 
belonging to one of two phases: 

• Static phase: if the grain is part of the sandpile, 

• Rolling phase: if the grain is not part of the sand- 
pile and rolls downward with a constant drift ve- 
locity. 



We consider the local slope 

Si = hi — hi+i 



(1) 



of the static grains as the variable controlling the dynam- 
ics of the rolling grains. Here, hi denotes the height of 
the sandpile at coordinate i, and we assume the pouring 
point at i = 1. 

At each time step a set of A^i small grains and N2 
large grains is deposited at the top of column 1 of the 
pile [see Fig. 1(a)] . These grains are considered to be- 
long to the rolling phase. In the case of thin flow we 
assume that the rolling phase is homogeneous and both 
types of species are mixed in the fluid phase and interact 
with the surface. Then one rolling grain of each species 
interacts with the surface at each time step, and can be 
converted from the rolling phase to the static phase. In 
the thick flow regime, the percolation effect is expected 
to take place, and the small grains screen the interaction 
of the large rolling grains with the pile surface. Thus in 
the thick flow regime, we consider that the large grains 
interact with the pile surface only when the number of 
small rolling grains at a given position falls below a given 
threshold e < A^2- 

The remaining rolling grains which do not interact 
with the pile surface are convected "downward" with unit 
"drift velocity" — i.e., they move to the next column at 
each time step. 

The dynamics of each rolling grain interacting with the 
sandpile surface is governed by its local angle of repose. 
We consider that the repose angle depends on the local 
composition on the surface, so we define Oap as the gen- 
eralized repose angle of a rolling grain of type a on a 
surface with local grains of type (3. We propose that 



to take into account the fact that large rolling grains roll 
more easily on top of small static grains than small grains 
roll on top of large static grains [since the surface "looks" 
smoother for large grains rolling on top of small grains, 
see Fig. ||(b)]. The repose angles of pure species Oaa lie 
between 621 and 612- 

Since a large grain roll easier than a small rolling grains 
on top of a layer of small static grains, we have 



hi < 



(3) 



Conversely, a large grain roll easier than a small grain on 
top of a layer of large static grains, so that 



722 < "12 ■ 



(4) 



The stratification experiments [p7[ use a mixture 
smaller less faceted grains and larger more faceted grains. 
Thus, the last inequality between the angles of repose 
is provided by the fact that the grains differ in shape 
of friction properties. Since the small species are the 
most rounded, then the repose angle of the smaller pure 
species is smaller than the repose angle of the large pure 
species — i.e.. 



'11 < "22- 



(5) 



Therefore, to mimic the experimental conditions for 
stratification |27|, we propose 



hi < 



< O22 < "12 ■ 



(6) 



We notice that the conditions (|^)-(^) are a consequence 
of the different size of the species, while the condition 
^11 7^ O22 is achieved for mixtures of grains with dif- 
ferent shapes, and does not depend on the size of the 
grains. Thus the model incorporates the size segregation 
and shape segregation in the definition of the angles of 
repose. In general the grains with the larger angle of re- 
pose will tend to be captured first. Thus the small grains 
will tend to be captured at the top of the pile, and the 
smoothest grains at the bottom of the pile. The perco- 
lation effects tend to segregate the small grains at the 
top of the pile, thus, it acts in the same way as the size 
segregation due to different angles of repose Eqs. (|)-(|). 



rolling 
grains 



721 < ^12 



(2) 
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FIG. 2. (a) Description of the discrete model, (b) The four 
dynamical angles of repose Oap depends on the composition 
of grains at the surface of the pile, and are chosen according 
of the four possible interactions between the two species of 
grains. As discussed in the text, the angles of repose satisfy 
O21 < Oil < 022 < 612, when species 1 are small and rounded 
and species 2 are large and rough, and 821 < 622 < dn < 6112, 
when species 1 are small and rough and species 2 are large 
and rounded. 

At each time step, the rolling grain interacting with 
the sandpile surface at coordinate i either will stop (by 
being converted into a static grain) if the local slope of 
the surface 

hi - hi+i < Sap = tan6'a/3, (7) 

where Sap is the local slope, or will continue to roll (to- 
gether with the remaining rolling grains) to column z + 1 
if 




FIG. 3. (a) Spontaneous Stratification: Results obtained 
with the discrete model for A''i = 10, = 10, e = 1, sn = 4, 
S12 = 7, S21 = 2, S22 = 5. Note the kink profile at which 
grains are stopped during an avalanche. The large rough 
grains are black and the small rounded grains are gray, (b) 
Segregation: Results obtained with the discrete model when 
the angle of repose of the large rounded grains is smaller than 
the angle of repose of the small rough grains for A'^i = 10, and 
N2 = 10, e = 1, sii = 5, si2 = 7, S21 = 2, S22 = 4. 



hi - hi+i > Sal3- (8) 

We iterate this algorithm to form a large sandpile of typ- 
ically 10^ grains. We assume that the substrate is made 
of a layer of large grains. 

The discrete model in the thick flow and thin flow 
regime gives rise to very similar results. Then in the fol- 
lowing we will concentrate only on the thick flow regime 
where most of the experiments are done, and where the 
grains are segregated in the rolling phase. Indeed, we will 
see that the differences between the thin flow regime and 
the thick flow regime arise only when cross-amplification 
processes are taken into account. Since the discrete 
model neglects these type of processes, then it gives rise 
to similar results for both flow regimes. 



B. Kink Mechanism 

Figure ^ shows the resulting morphology predicted 
by the discrete model when the large grains are rougher 
than the small grains. This fact is quantified by taking 
into account that the angle of repose of the pure large 
species is larger than the angle of repose of the pure small 
species, 622 > On. The stratification is qualitatively the 
same as that found experimentally [^7|] , both in regard 
to the statics of the sandpile [seen in Fig. ||(a)], and also 
in regard to the dynamics (seen in Fig. Q). 

The mechanism for spontaneous stratification involves 
the formation of a kink at which the grains are stopped 
during an avalanche. After a pair of static layers is 
formed with a layer of large grains on top of a layer 
of small grains, the angle of the sandpile is close to 6122 
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[Fig. 0(a)]. Since the surface of the sandpile is made of 
large grains and 6*22 < 0i2, a thin layer of small grains is 
trapped on top of the layer of large grains as more grains 
are poured onto the cell. These small grains smooth the 
surface without changing significantly the sandpile an- 
gle, and allow rolling small grains to go further down 
(since On < 6*22) [Fig- §(b)]. The large grains are roUing 
down on this thin layer of small grains without being 
captured (^21 < ^22), and are the first to reach the bot- 
tom of the sandpile, giving rise to segregation [Fig. 0(b)] . 
When the flow reaches the base of the pile, the grains de- 
velop a profile which displays a kink where the grains are 
stopped: the small grains are stopped first at the kink 
since 6*21 < 6*11, and a pair of layers begins to form, with 
the small grains underneath the large grains [Fig. 0(c)]. 
The kink moves upward (in the opposite direction to the 
flow of grains) until it reaches the top of the pile and a 
complete pair of layers has been formed [Fig. 0(d)] . 




If, on the other hand, we consider ^22 < ^11 in the 
model — corresponding to a mixture of smaller "more 
faceted" grains, and larger "less faceted" grains — we find 
the complete segregation of the mixture but no strati- 
fication [Fig. 0(b)]. In this case, the small and faceted 
grains are to be found near the top of the pile, while the 
large and rounded grains are found near the bottom of 
the pile. Thus, the control parameter for spontaneous 
stratification appears to be the difference in the repose 
angle of the pure species, which quantifies the difference 
in shape of the grains. 



As seen in Fig. 0(a), before the layers appear there 
is an initial regime where only segregation is found. At 
the onset of the instability leading to stratification, a few 
large grains are captured on top of the region of small 
grains near the center of the pile where the angle of 
the pile is 9 ~ 6*11. The repose angle for large rolling 
grains is now ^22- Thus, if ~ On < 622, more large 
grains can be trapped (since the angle of the surface is 
smaller than the repose angle), leading to the first sub- 
layer of large grains and then to stratification. On the 
other hand, if 2± On > 6'22, no more large grains can get 
captured, the fluctuation disappears, and the segregation 
profile remains stable. Thus this picture suggests that, 
in agreement with [^, the segregation profile observed 
in the initial regime is "stable" so long as ^22 < ^11, and 
unstable (evolving to stratification for large enough sys- 
tems) when Oil < 622- This explains qualitatively why 
the control parameter for the stratification-segregation 
transition is the angle of repose of the pure species. 



FIG. 4. Dynamics obtained with the discrete model, show- 
ing the formation of a pair of layers at the kink. Notice the 
kink moving uphill at constant velocity U|. 



III. CONTINUUM THEORY FOR GRANULAR 
FLOW 

A theoretical approach to segregation in granular fiow 
is not a simple task. On one hand, extensive numerical 
simulations (such as molecular dynamics [0,0,^-|3|, or 
lattice gas models p^-p7[]) have been able to reproduce 
several segregation phenomena observed experimentally. 
On the other hand, an analytical approach of segrega- 
tion, considering granular medium as continuum, would 
be instructive, since it would allow to reduce the problem 
to a small set of parameters that control the system, and 
then to get a clear phenomenological understanding of 
the problem. However, an important question is whether 
the grains can be treated as a continuum medium. A 
continuum approach means that we will be able to re- 
place the grains belonging to a same region of space by 
some average quantities, for example the density, or the 
averaged speed. But specific grains can play a peculiar 
role, i.e., in the case of arching effects, where an arrange- 
ment of several grains in a special way (an arch) supports 
the above grains, prevents the flow by gravity, and cre- 
ates large discontinuities for the local density or the force 

In this section, we study spontaneous stratification us- 
ing a continuum formalism. We will show that the results 
obtained are very close to the experiment, and then we 
extract the important parameters of the problem. Our 
results shows that a continuum formalism is able to re- 
produce a complex phenomenon in granular matter such 
as spontaneous stratification. 

The continuum theory of granular fiow for a single- 
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species sandpile was proposed by Mehta et al. and 
Bouchaud et al. |34|. Bouchaud et al. recognized the 
necessity of treating differently the rolhng grains (rolhng 
on the surface of the surface), and the static grains that 
form the sandpile. The rolhng grains are in a "liquid" 
state, flowing down by gravity and interacting with the 
static grains in the "solid" state. The dynamics of the 
sandpile is then given by how the rolling grains move on 
the surface and how they interact with the sandpile. The 
equations for those two quantities for the case of single- 
species sandpiles are 



dt 



R[x{t)] = r(i?, ft). 



dh{x, t) 
dt 



-T{R, h). 



(9a) 
(9b) 



Here R[x{t)] is the height of the layer of rolling grains 
in the Lagrange representation (following the rolling 
grains), and h{x,t) the height of the sandpile [Fig. |l|(a)]. 
The rolling grains interact with the sandpile through 
r(i?, /i), which represents the height of rolling grains be- 
coming part of the sandpile. Equations (^ conserve the 
number of grains. 

We need to define how the rolling grains move on the 
surface, and how they interact with the sandpile. Two 
forces are acting on the rolling grain: (i) gravity and (ii) 
friction coming from the collisions of the rolling grain 
with the grains of the sandpile. The number of collisions 
at the scale of interest for continuum equations (several 
times the size of a grain) is large enough to say that at 
this scale, the speed is always at its limit value where the 
two forces balance exactly, and the motion of a grain is 
overdamped. Equations ( pa) ) and (9b) can be rewritten 
in Eulcr representation (looking at the height of grains 
at a given position x) 



dR[x,t) 

di 
dh{x, t) 

dt 



dR{x,t) 
dx 

-r{R,h). 



r(i?, h), 



(lOa) 
(10b) 



Here R{x,t) is now a function of the time and the posi- 
tion X, and we assume the pouring point at a; = 0. The 
rolling grains are moving down at a speed v{x,t) in the 
positive X direction, where v{x,t) could depend on the 
local slope of the sandpile. However, in the spirit of cap- 
turing only the essential physical mechanisms to recover 
the observed phenomena, we will consider the speed v to 
be constant all over the sandpile. 

Bouchaud et al. proposed a form for r(i?, h) in the 
general case, that has been simplified by de Gennes |3^] 
for the simpler case of a continuous flow of rolling grains. 
In this case, there is no discrete avalanching and only the 
repose angle must be included in the formalism | |5l| ] . The 
interaction term r(i?, h) is 



T{R,h) ^-f[9{x,t) -er]R, 



(11) 



d{x,t) 



dh{x, t) 
dx 



(12) 



is the local angle of the surface (we will make no distinc- 
tion between the angle of the surface and the tangent of 
the angle). 

Equation ( pi] ) states that the rate of the interaction is 
proportional to the number of grains interacting with the 
sandpile. Rolling grains will become part of the sandpile 
if the angle of the surface 9{x, t) is smaller than the re- 
pose angle 9r ("capture"), while static grains will become 
rolling grains if9{x, t) is larger than 9r ( "amplification" ) . 
The constant 7 is the inverse of the time on which this 
interaction is effective because the distance on which a 
rolling grain is stopped when 9{x,t) is smaller than 9r is 
v/"f. This distance must scale with the size of the grain 
d. Thus EEA 



u/7 ~ d. 



(13) 



The linear dependence of T{R, h) on the quantity 
[9{x,t) — 9r] can be understood as a first order devel- 
opment of a more complicated function of 9{x,t). As 
soon as the angle of the sandpile is far from the repose 
angle, non-linear terms must be added (avoiding the non- 
physical divergences that could be found with the linear 
development). Finally the proportionality of r{R,h) to 
i?(x, t) is meaningful in the case where all the grains in- 
teract at each time with the surface (thin-flow limit). 
This imposes the height of the rolling layer to be of the 
order of the grain size R{x,t) ~ d. However we will ar- 
gue that this approximation might be still valid in the 
case of thick flows as well, since the interaction might be 
proportional to the pressure exerted by the fluid phase, 
which in turn is proportional to R{x,t) for a fluid p^ . 

To treat the problem of segregation in granular flow 
composed of binary mixtures, Boutreux and de Gennes 
(BdG) have extended this formalism to the case of 
flows made of two types of grains, identifying three vari- 
ables: the two heights of rolling grains Ra{x,t) (i.e. the 
total thickness of the rolling phase multiplied by the lo- 
cal volume fraction of the a grains in the rolling phase at 
position x), and the height of the sandpile h{x,t). BdG 
generalized Eqs. (fol) to 



dRa 

dh 
'di 



dRc , ^ 
dx 



a 



1,2 



(14a) 
(14b) 



where the interaction term Fq takes into account the 
conversion of the rolling grains into static grains, and is 
deflned through the collision matrix M^^ 



where 



13=1 



(15) 



7 



The colhsion matrix, governing capture and amplifica- 
tion, depends on the local angle of the pile 9{x,t) and 
on the local composition of the surface of the sandpile 
4>a{x, h), which is a function of x and h. However, writing 
an equation of evolution for (j)a{x, h) is not easy. When 
both rolling species are captured, the height of the sand- 
pile increases, and 4>a{x, h) is given by 



the thin and thick flow regimes, in order to understand 
stratification and segregation as seen in the experiments. 



dh/dt' 



and 



(16) 



(17) 



As soon as amplification dominates for one or both 
species, Eq. ( |l6|) is no longer valid. When Fq > for 
both species, the height of the sandpile decreases, and 
(j)a{x,h) does not have to be updated. Finally, when 
grains of type a are captured and grains of type f3 are 
amplified, the composition will be 4>aix, h) = 1. 

The general form of the collision matrix is defined by 
taking into account a set of binary collisions between a 
rolling and a static grains 



^.r ^ ( ai{0)(t)i ~ bi{e) X2\u)ipi 

' XI (0)02 a2{e)4>2-h2{e) 



(18) 



This definition involves a set of a priori unknown colli- 
sion functions contributing to the rate processes: aa{9) 
is the contribution due to an amplification process (when 
a static grain of type a is converted into a rolling grain 
due to a coUision by a rolling grain of type a), ba{9) 
is the contribution due to capture of a rolling grain of 
type a (when a rolling grain of type a is converted into 
a static grain), and Xa{0) is the contribution due to a 
cross- amplification process, (the amplification of a static 
grain of type /3 due to a collision by a rolling grain of 
type a). 

BdG proposed an analytically tractable form for 
the collision matrix which includes capture and amplifi- 
cation, in the case where the two species have the same 
size but differ in respect to their repose angle. BdG found 
the steady state solution in the geometry of the silo, in 
the case where the repose angle does not depend on the 
composition of the surface. The BdG model explains 
the segregation of the two type of grains in different re- 
gions of the sandpile, but not the phenomenon of spon- 
taneous stratification. A generalization of the model to 
include also the different sizes of the grains shows a 
smooth segregation of the species with the concentrations 
of static grains behaving as a power law of the position 
along the surface pile; results valid only when the species 
do not differ too much (a fact that allows to perform lin- 
ear approximations of the collision functions around the 
angles of repose ) . In the following we propose a form 
of the collision functions valid when there is a large differ- 
ence in size and shape between the species and we treat 





(b) 



segregation 



FIG. 5. Dependence of the repose angle for the two types 
of rolling grains on the concentration of the surface of large 
grains 4>2. (a) An essential ingredient to obtain segregation is 
that 6*22 > On, (b) while for segregation we require > 6122. 
For the numerical integration, we use the linear interpolation 
between 02 = and 02 = 1 as plotted here. 



IV. CONTINUUM MODEL FOR 
STRATIFICATION OF GRANULAR MIXTURES 

A. The collision matrix and the generalized angle of 
repose 

We present the collision matrix Map that includes the 
effects of capture and amplification, and the dependence 
of the repose angle on the composition of the surface as 
discussed in |37|. Considering one type of rolling grain a, 
we define capture and amplification as follows: if the lo- 
cal angle of the surface 9{x,t) — —dh{x,t)/dx is smaller 
than the generalized angle of repose of the rolling specie 
a da , the a rolling grains will be captured. In the reverse 
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case 6{x,t) > 9a, the static grains at position x will be 
amplified. Amplification will cause locally a small hole, 
i.e., the two types of static grains will be amplified the 
same way according to their local concentration 

The generalized repose angle 9a of each type of rolling 
grain is now a continuous function of the composition of 
the surface 9a — 9a{(t>p) (see Fig. ||) 



92U2) 



m(p2 
m02 



711 

921 



-m0i + 022, 



(19) 



where m = 612 — 9ii = 622 — ^21- We have assumed the 
difference 



^ = 9i{cj)2) ~ 02{h) 



(20) 



to be independent of the concentration (j)2. For simplic- 
ity, we choose a linear interpolation between the extreme 
values da{(f>i3 = 0) and 9a{4)p = !)• If we label the small 
grains with a — 1 and the large grains with a — 2 
and denote the extreme values of the repose angle by 
0a{4>l3 = 1) = ^ap as in Sec. ||, the difference in size 
implies [see Fig. §(b)] 



712 < ^21. 



(21) 



If species 1 are the smallest then 9i{(p2) > 6*2 ('^'2) for 
any composition of the surface (j)2 — i.e., the small grains 
are always the first to be captured. Thus, the angular 
difference ip is determined mainly by the difference in 
size between the species, and the value of ip determines 
the degree of segregation; the larger tp, the stronger the 
segregation. The repose angles of the pure species 9ii 
and 622 will be intermediate between 612 and 622, and 
their relative values will depend on the shape of the two 
species as discussed in Sec. [I A . 

For mixtures of grains with different shapes or fric- 
tion coefficients we have 9ii ^ 622, and 9i2 — 621 if the 
species have the same size. If 6*11 is the repose angle of 
the pure round species and ^22 is the repose angle of the 
pure rough species, then < ^22- If the species have 
the same shape or friction coefficients then 9ii = 622- 

When the size ratio is close to one (^2/^1 ~ 1-5), the 
angle ip is expected to be small and then it is plausible 
to linearize the collision functions around the angles of 
repose as done in . When the is a larger difference in 

size {d2/di ^ 1-5), we expect ip to be large, and we ap- 
proximate the collision functions and define the following 
collision matrix: 



Ma0 = 



where 



711,501 e,,„[50i, cPi] 7i2'502e„_[^^2, M 

721^eie„_[501, 02] 1226920^392, <t>2] 



9a=9{x,t)-< 



(22) 



(23) 



Here, jaa > has dimension of inverse time, and v/jaa 
represents the length scale at which a rolling grain will 
interact significantly with a surface at an angle above or 
below the angle of repose. 

The functions Q{S9a, 4>a) distinguish capture from am- 
plification. Ggeit treats the case of self amplification, i.e., 
amplification of static grains a by rolling grains a, and 
©cioaa the case of cross-amplification, i.e., amplification of 
static grains a by rolling grains (3. More specifically, we 
consider 



K) = 



1 



if 9-9,, 
if 9-9„ 



< 
> 



(24) 



These equations mean that in the case of capture, the 
amount of rolling grains of type a converted to static 
grains is proportional to the number of grains interact- 
ing with the surface, i.e., Ra', in the case of amplification, 
the amount of static grains converted to rolling grains is 
proportional to Ra again and to the number of static 
grains a on the surface, i.e., (pa- Accordingly, we obtain 



e.rU0-Oa,<Pp) 







if 
if 



< 
> 



In terms of the collision functions defined in (|TE 
have 



(25) 



we 



aa{0) = 7aa U[9{x,t) - 9a{^p)], 
ba{0) = lac n[0„(00) -0(a;,t)], 
X(i{9) = 70„ Xl\9{x,t)-9p{<Pp% 



where 



if a; < 
if a; > 



(26a) 



(26b) 



This form for the matrix Map supposes that amplifi- 
cation (respectively capture) occurs only when the local 
angle is larger (respectively smaller) than the repose an- 
gle. Thus, for a given angle only capture or amplification 
can occur at a given time; an assumption which is ex- 
pected to be valid when the grains differ appreciable (^/^ 
large). This fact gives rise to strong segregation effects 
with exponential behavior of the concentrations as we 
will corroborate when calculating the steady state solu- 
tion of the problem in Sec. V A. In contrast, the form of 
the collision functions proposed in |36 53] assumes that 



amplification and capture are linear functions of the lo- 
cal angle. For a given angle both amplification and cap- 
ture act simultaneously, and the repose angle corresponds 
to the angle where amplification and capture are equal. 
This linear approximation is valid when the species do 
not differ much (for small V'), and gives rise to a weaker 
segregation of the species — with a power law behavior of 
the concentrations as a function of the position along the 
pile surface — than the segregation found in this study 
which is valid when the species differ appreciable in size 
and shape. 
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We focus here on the dependence of the repose angles 
on the composition of the surface. We will consider that 
the other parameters — i.e., ^acx and Va — do not depend 
on the composition of the surface and on the species con- 
sidered 



7ii — 722 = 7j 



(27) 



In the conditions of the experiment [p7| , where an equal 
volume of the two species is poured at the left side of the 
cell, the boundary conditions are 



i?«(0,i) = i?° = 



(28) 



These equations are meaningful if all the rolling grains 
interact each time with the surface, and do not interact 
directly with each other. Then, when there is perco- 
lation in the rolling phase, the interaction terms have 
to be modified accordingly. In the next section we will 
show how to introduce the percolation effect characteris- 
tic of the thick flow regime. However, we will show that 
the above equations — although strictly valid in the thin- 
flow limit, i.e., Ra — da — will be also needed to describe 
completely the thick flow regime. 



B. Thick flow regime and the percolation effect 

An important consideration of the model proposed 
above is that the two types of rolling grains do not in- 
teract significantly each other when the flow of grains is 
small, i.e., Ra(x,t) < e, where e is a cut-off of the order 
of the grain size da- This hypothesis may not be valid 
in experiments where the input flux is large, where the 
layer of rolling grains could be larger than one grain size. 

In the case of a thick layer of rolling grains, the phe- 
nomenon of percolation occurs; during the flow of grains 
down an inclined plane, the small rolling grains percolate 
downward through the gaps left by the motion of larger 
rolling grains |2^j2^,^ [Fig. |l|(b)]. Therefore, as they are 
convected down, the rolling grains segregate, the concen- 
tration of small grains being large at the bottom of the 
rolling layer, and the concentration of large grains large 
at the top of the rolling layer. As only the lowest grains 
of the rolling layer interact with the sandpile, the above 
effect imposes that small grains will preferentially inter- 
act with the static grains. We include the percolation 
effect in the continuum model and see how the results 
related to spontaneous stratification are affected. 

To model this effect, we consider the extreme situation 
where the segregation is present as soon as the grains 
start to roll, i.e., up to the top of the sandpile, and that 
only the small grains will interact with the sandpile as 
soon as Ri{x,t) > e, with e <C e is thus the height 
of small grains below which the percolation stops to be 
effective, and for which the large rolling grains start to in- 
teract with the surface. The interaction terms discussed 




in Sec. IV A have to be modified accordingly. 



FIG. 6. Thick flow regime: Morphology resulting from 
the numerical integration of the continuum equations when 
the percolation effect is present, (a) Stratification of a mix- 
ture of large rough grains and small smooth grains. The pa- 
rameters used are RI — R2 — 1, e = 0.25, tan(Sii) = 0.5, 
tan(6'22) = 0.6, tan(6'i2) = 0.9, tan(02i) = 0.2, 711 = 722 = 1, 
721=712 = 0.1, and vi = V2 = 1. The dynamics obtained by 
the numerical integration of the equation of motion proposed 
in the text show the formation of a pair of layers through 
the kink mechanism. We also see that the rolling grains are 
stopped at the kink in similar fashion as was observed in the 
stratification experiment |2^ and in the discrete model in Fig. 
^ (b) Complete segregation of a mixture of small rough grains 
and large smooth grains obtained with the same parameters 
as in (a) except for the angle of repose of the pure species 
tan(Sii) = 0.7, tan(^22) = 0.5. Notice the total segregation 
of the mixture. The region of mixing is concentrated in a 
small region in the center of the pile. 



This model implies that except for a small region where 
Ri < e, the surface of the sandpile is always divided in 
two regions. The upper region corresponds to the part 
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of the surface where small roUing grains are present. As 
they are the only one to interact with the sandpile, the 
sandpile will be made only of small grains. In the lower 
part, only large grains are present (except for the e part 
of small grains). The division point between the two re- 
gions is able to move with time. 

For the upper part of the sandpile, defined by 
Ri{x,t) > e, only small grains are captured, with 
R2ix,t) = Rl = i?°/2 and (l)i{x,t) = 1. The equations 
become 



dt 
dh 



dx 



_ ^ -^^0 ^ OnKRi + R2) 



(29a) 
(29b) 



When Ri{x,t) < e, the percolation effect disappears 
and the equations of evolution of Ra and h can be con- 
sidered to be the same as the one defined by the collision 
matrix ( |2^ ) valid for the thin flow limit. 

The effectiveness of the interaction (given by 7 before) 
must be multiplied now by the vertical pressure acting 
on the lowest rolling grain interacting with the surface, 
i.e., the weight of the column of rolling grains above the 
interacting grain |p^ . This explains the presence of the 
R2{x,t) terms in the interaction term Fi in Eq. (p^. 

We numerically solve Eqs. dl4) with the interaction 
terms given by ( |2^ ) when Ri{x^t) < e and by ( [29| ) when 
Ri{x,t) > e as discussed above. We find qualitatively 
the same results as in the experiments [ pTf , with stratifi- 
cation appearing when the large grains are rougher than 
the small grains [622 > ^11, Fig- ||(a)]. The dynamics 
of stratification show also the formation of the kink as 
observed in the experiments and in the discrete model 
shown in Fig. ^. When the small grains are rougher than 
the large grains (^n > 022), then we obtain the strong 
segregation of the mixture with the small rough grains 
found at the top of the heap (Fig. ||(b)) and a small re- 
gion of mixing in the center of the pile |53] . 



C. Thin flow regime 



we find that the transition to segregation does not oc- 
cur sharply at (5 = but occurs at a small negative 
value Sc < 0, which depends on the value of the cross- 
amplification rates 70/3. When the cross-rates are zero, 
then the transition occurs at ^ = as we found for the 
thick flow regime above. For smaller values oi 6 < Sc < 
we find the complete segregation pattern found in the 
experiments. Thus, the presence of cross- amplification 
processes — which appear only in the low flux limit — 
shifts the transition from stratification to segregation to a 
value Sc different from zero |Q. However, the stratifica- 
tion we find for (5c < (5 < is less pronounced than in the 
case 5 > 0: the layers do not go to the top of the sand- 
pile, and the layer of small grains is very thin. Similar re- 
sults at the neighborhood of the stratification-segregation 
transition have been found with a microscopic model of 
grains dynamics [^5| . We also find a kink, correspond- 
ing to the growth of the new pair of static layers, with a 
well-defined steady-state profile and upward velocity for 
S>0. 

This model, valid for thin flows, is qualitatively close 
the preceding one valid for thick flows, where we also 
found a regime of complete segregation separating the 
surface of the sandpile in two regions. This can be un- 
derstood noting that the dependence of the repose angle 
on the composition of the surface in the thin flow model 
is comparable to the effect of percolation in the rolling 
layer in the thick flow regime. Since for a given compo- 
sition (f)p, the repose angle of the small grains is always 
larger than the repose angle of the large grains (Fig. ||), 
the small grains are also the first to be trapped before 
the large ones. The composition of small grains on the 
surface increases, amplifying the effect, up to the point 
where the surface is made of only small grains. 



V. ANALYTICAL RESULTS 

The thick flow and the thin flow regime show similar 
results as long as there is a large difference in size be- 
tween the grains. We will now show analytical results 
valid for both regimes. 



In the absence of the percolation effect, Eq. ( |2^ ) re- 
mains valid along all the pile since in the small flow limit 
a thin rolling phase is expected with Ra{x,t) < e. How- 
ever, we assume a large difference in size, so that strong 
segregation effects are expected anyway, and the collision 
functions are expected to behave as in (^2|). We then 
solve Eqs. (p^)-(p2[) numerically. Defining the "control 
parameter" 

<5 = 022-^ii, (30) 

we obtain stratification when 6 > 0, i.e., when the large 
grains are rougher than the small grains as in the ex- 
periments and in the thick flow regime above. However, 



A. Steady state solution 

We next study the stationary solution of Eqs. ( p^ ) and 
then we study how a perturbation to this solution can 
be amplified giving rise to stratification. For simplicity, 
we consider the geometry of a silo of lateral size L, i.e., 
all the rolling grains are stopped when they reach a wall 
at position x — L [Fig. [|(a)]. Moreover, we assume the 
difference 

^ = e2{<t>2) - Bl{4>2) (31) 

to be independent of the concentration 02 • We seek for 
a solution where the profiles of the sandpile and of the 
rolling grains are conserved in time. Thus, stratification 
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that corresponds to periodic variations in time of the dif- 
ferent variables of interest cannot be observed for this 



solution. The conservation of the grains gives [56 

V R° 



and we impose 



dh 
'dt 



L 



(x) = 0. 



(32) 



(33) 



We assume that, as in the experiment, an equal volume 
mixture is used, so 



and the boundary conditions are 



RciL) 

i?a(0) = 



= 0, 
i?72. 



(34) 



(35) 



The steady-state solution of ([l^)-(^) shows total seg- 
regation: except in a region of size of v at the center 
of the sandpile where the grains are mixed, the sandpile 
is exclusively made of small grains in the upper part of 
the pile, and of large grains in the lower part of the pile. 
The details of the calculations can be found in Q . Here 
we present a simplified results for the case of v/j <C L, 
since we expect (p^. 

At the upper part of the pile, for < a; < Xm, 
with Xm = L/2 — v/l^fip), only small grains are present 
{4>i{x) — 1, and (j)2{x) = 0), and the profiles are 



R2{x) = i?72, 
(x) - 6*11 - 



i(l+72i/7)/2-x' 



(36a) 
(36b) 
(36c) 



At the lower part of the pile {xm < a; < L), we find that 
after a small region of size of the order of v/{'jtp) (or 
equivalently, of the order of the size of the grain), mainly 
large grains are present. The profiles are 



01 (x) 
Ri{x) 



exp 
2v 



9{x) - 022 = -m<f)i{x) - 



\X Xm } 

V 

(l)i{x)R{x), 



l{L - x)' 



Here 



R{x) = Ri{x) + R2{x) = i?°(l - x/L), 



(37a) 
(37b) 
(37c) 

(38) 



and m = 022 - 02i = (^12 - ^ii- 

When the system is in this stationary state, the an- 
gle of the upper part of the sandpile (where only small 



grains exist) is almost On; without cross-amplification 
(721 = 0), the sandpile would be at its repose angle 0ii, 
and because of cross-amplification, the angle is signifi- 
cantly reduced. The divergence around x = L/2 insures 
that all the small grains are captured despite cross am- 
plification. Moreover the angle for the low part of the 
sandpile is almost 022- If we now suppose a fluctua- 
tion in Ri such that few small rolling grains enter the 
low part of the sandpile, capture will be of the order of 
7(6*22 — 6i{(l>2))Ri, and amplification will be of the or- 
der of 7(^22 — d2{<f>2))R2- As i?i <C i?2, amphfication 
will dominate capture, and the small rolling grains will 
go downhill up to the point when capture starts to ex- 
ceed amplification. As soon as this small layer of rolling 
grains exist, the large grains are convected down with- 
out being capture. This excess of large grain creates a 
kink at the bottom of the pile going uphill where large 
grains are now stopped; the first layers are formed and 
the same process starts again. This may explain why 
when cross-amplification terms are present the transi- 
tion stratification-segregation is shifted from 6 — to 

5^Sr<0. 




I' 



stable 
' segregation 



stratification 




022>6i 



FIG. 7. Stability analysis of the steady state solution of the 
equations of motion for granular flow of mixtures in the ge- 
ometry of the silo. When 5 < we find that the steady-state 
solution is stable, so that this solution is the physical solution 
when S < 0. On the other hand, when 5 > the steady state 
solution is unstable and it evolves into stratification as in the 
experiment. 
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When percolation is present this shift disappears since 
cross- interactions are greatly reduced by the screening ef- 
fect of the small rolling grains at the bottom of the rolling 
phase. 

To analyze the stability of the steady-state solution for 
the different phenomenologic parameters, we impose the 
steady-state solution as the initial condition, and then 
we look numerically for the stability of the profile under 
perturbations (Fig. ^). For On > 6*22, the steady state 
solution is stable, so that Eqs. ( ^6| ) and ( p7| ) are the 
physical solution for this case. In this case, only segrega- 
tion is observed, and the sandpile conserves in time the 
profiles (|3^ ) and {^t\). For 9ii < 622, the steady-state 
solution is unstable (evolving to stratification), just as in 
p7| . The onset of the instability is clearly seen, where 
the kink first start at the center of the pile where some 
large grains are captured on top of small grains as seen 
in Fig. for 6*22 > dn. 

Our model predicts the complete segregation of the 
species when 622 < dii- The concentration of small 
grains drops very fast near the center of the pile in a 
characteristic region of the order of v/{jip) ~cm accord- 
ing to experiments As discussed in Sec. IV A this 
is a consequence of our assumption of large difference in 
properties of the species. 



B. Analytical Shape of the Kink 



du 



dfiu) 



du 




(40) 



with 



^11 = 



'11 



^21 — ^21 ^ I 



(41) 



We obtain the shape of the low part of the kink. For 
u < 0, f{u) = and for m > 0, f{u) and Ri{u) obey the 
following equations: 



9/ ^ ('?ii+^2i)i^72-^;(5n/(^) 
du V]-/^ + wf{u) — 

Ri{u) = -wf{u) + R^/2 



(42) 



where w = / {v + v-\). Then the lower part of the kink 
is characterized by a linear dependence 



f{u) (X 



(^ll+<52l)/2 



(43) 



Finally, this solution is valid up to the point where the 
large grains start to be captured, i.e., df /du = —52i- 
From Eqs. (^) and (|3|) we obtain the following inequal- 
ities: 



We saw in a preceding section that the layers were con- 
structed through a kink where all the grains are stopped. 
A precise analysis of this kink would allow us to under- 
stand the dynamic of the process and the characteristic 
size of the layers. 

From the simulations we observe that the kink has a 
well define uphill speed and that its shape is preserved 
in time. To obtain an analytically tractable set of equa- 
tions, we then look for a possible steady-state solution for 
the shape of the kink, and thus for R{x,t) and h{x,t). 
We also must make the assumptions that far below and 
above the kink, the sandpile has a constant angle a-nd 
is only made of large grains. These assumptions are ver- 
ified by comparing the obtained result for the shape of 
the kink with the numerical results. 

The existence of a stationary solution for the kink im- 
plies that Ra{x,t) and f{x,t) = h{x,t) 4- Oqx are func- 
tions only of 



X + v^t, 



(39) 



where U| is the uphill speed of the kink. For the lowest 
layer of the kink composed of small grains (the upper part 
of the pile), as only small grains are captured {(j)i{u) — 1, 
R2{u) = R°/2), Eqs. ( pl| ) reduce to equations for Ri{u) 
and f{u) (we assume 721 = 7), 



% > 



<7i?° 



(44) 



The solution of the equations for the upper layer of the 
kink (the lowest part of the pile) where only large grains 
are present can be obtained in the same way and is 



(45) 





1 / Wj622u\ 











where 



^2 



(46) 



We then find that the shape of the upper part of the 
kink is exponential, and that the stationary solution ex- 
ists only for S22 < 0. 

Thus we see that the existence of the stationary solu- 
tion for the kink implies that 



711 



721 



< 



(47) 



and that the sandpile is built on an intermediate angle 
between those two extreme values. We are now able to 
distinguish two types of stratification. When Eq. (^^ is 
satisfied, a stationary solution exists for the kink and the 
layers are well formed. When Eq. (47) is not satisfied, 
the layers are very asymmetric, with thin layers of small 
grains. 
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C. Wavelength of the Layers 



The layer thickness A is defined by the width of the 
kink, A = lim„^oo /("). From Eq. (^Sf), we obtain 
A = /w, which is a consequence of the conservation 
law stating that all the rolling grains are stopped at the 
kink. This relation is obtained assuming that the density 
of the fluid phase is the same as the density of the bulk. 
However, in general we have ^n-^id < Mbuiki where /ifluid 
and /ibuik are the number of rolling grains per unit vol- 
ume of the fluid phase and bulk, respectively. Then we 
obtain 



^ ^ Mfluid {v + V]) ^0 

A*bulk 

Furthermore, Eq. (^3|) implies that 



(48) 



(49) 



where Cf is a numerical constant that does not depend 
on w, 7, or i?". Then we obtain 



A- 



thin 



fluid 
Aibulk 



V 



R° 



[thin flow regime] 



(50) 



This relation is relevant to the thin flow regime where 
w/7 ~ d and i?° ~ d are both of the order of the size of 
the grains, and the velocity of the rolling grains is con- 
stant since the rolling phase is homogeneous. However in 
the thick flow regime the mean value of the velocity of 
the grains scale with the thickness of the rolling phase as 
found recently experimentally [Eq] 



(51) 



where C is a numerical constant that does not depend 
on 7 or i?'^, but may depend on the angles of repose and 
other features of the grains. In this case Eq. ( |4^ ) be- 
comes 



At. hi 



tliick 



A*fluid 
Aibulk 



C 

1 H I i?" [thick flow regime] . 



(52) 

Thus we see that we expect a linear dependence of the 
wavelength of the layers as a function of the thickness of 
the rolling phase, i.e. as a function of the flux of grains 
(if the plate separation is maintained constant |28| ) , but 
the proportionality constant is different in the thin flow 
regime and the thick flow regime. From ( ^ ) and ( |5^ ) 
we can obtain an approximate estimation of the value of 
the cross-over from the thin flow regime to the thick flow 
regime i?° assuming that Xthin{Rc) ~ Athick(^c); then 
Rc ~ I'thin/ (C7), where ?;thin is the velocity of the rolhng 
grains in the thin flow regime. 

Typical experimental values of the phenomenological 
constants of the problem for a system of i = 30 cm. 



di = 0.27 mm, ^2 = 0.8 mm, are u w 10 cm/sec, 
7 « 20/sec, tan?/; = tan 6*11 -tan 6*21 ~ 0.2, /ifluid/z^buik = 
0.85, and C/Ct « 1. Therefore, v/{"ftainp) ~ 2.5 cm. 
We also obtain « 0.5 cm as the value of the cross-over 
flux from the thin flow regime to the percolation regime. 
The above predictions regarding the dependence of the 
wavelength on the flux of grains have been confirmed in 
recent experiments [Eq,E3]. 



VI. DISCUSSION 

In summary, we develop a mechanism to explain the 
spontaneous stratification observed in j2^. This mecha- 
nism is related to the dependence of the local repose angle 
on the local surface composition. When we consider only 
capture (corresponding to the cellular automaton model), 
we find that spontaneous stratification occurs only when 
the repose angle of the large grains is larger than the re- 
pose angle of the small grains {622 > corresponding 
to large grains rougher than small grains). This result is 
in agreement with the experiments. Stratification is also 
obtained when 6*22 is slightly smaller than On as soon 
as we include cross-amplification process, as we find in 
the thin flow regime with the continuum model. When 
O22 < On, the model predicts almost complete segrega- 
tion, but not stratiflcation. These results are in agree- 
ment with experiments ]2^ ]. 

It is interesting to compare those results with a pre- 
vious study 1 27 1 where a continuum model was used ne- 
glecting cross- amplification (i.e., amplification of a static 
grains by P rolling grains). In this case, it was shown 
that amplification never occurs during the formation of 
the sandpile, and stratification is obtained in the same 
conditions as in the experiments. 

We discuss the possible importance of percolation ef- 
fects in the rolling layer, that takes place when the input 
flux of grains is large. Although the continuum model 
must be deeply modified, the results are surprisingly close 
to the thin flow model. 

The model describes well the static picture of the sand- 
pile of with alternating layers made of small and large 
grains, and also reproduces the dynamics, where the lay- 
ers are built through a kink mechanism. The numerical 
simulations suggest that the motion of the kink is sta- 
tionary. It allows us to make quantitative predictions for 
the dependence of the size of the layers on the different 
parameters of the problem. 

From the theoretical point of view, it is the flrst time 
that this continuum formalism for granular flow is di- 
rectly compared to experiments. The quality of the re- 
sults suggest that this formalism includes the essential 
features of the physics of granular flow. Using this for- 
malism, it may be possible to have theoretically access 
to 3D problems, such as periodic segregation in rotat- 
ing cylinders p6[ , that are out of reach with present-day 
computer simulations. 
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